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АВ5ТКАСТ 


Ісі бпесіе describes four microcomputer programs for 
evaluating sonobuoy effectiveness. The programs are based 
the Gaussian Signal excess model. The first program 
can be used to compare different sonobuoy employments, for 
example, a comparison between deep and shallow hydrophone 
depths. The second program gives a user the opportunity 
ШО display the effective extent of his pattern. The third 
program provides localization information using a buoy 
pattern and contact status. The program is based on a 
model that uses both positive and negative information 
from a field of sonobuoys to estimate target position; 
and is designed to give the operator an indication of where 
to deploy additional sonobuoys in order to convert single 
buoy to multiple buoy contact. Topographic mapping 
procedures are used to display the information in a graphical 
format that could prove useful to antisubmarine warfare 


aircrews. 
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I. BACKGROUND 


The localization problem facing an antisubmarine warfare 
(ASW) air crew is very complex. The crew must properly 
employ a variety of acoustic and non-acoustic senso-s and 
accurately process the resultant information tc properly 
classify all targets and to identify, localize and pros=cute 
mae target of interest. The purpose of the research that 
1s reported in this thesis is to develop techniques to 
analyze the effects that a changing environment can have 
on this activity. The techniques that were developed are 
relevant to planning sonobuoy patterns for specific 
scenarios and they also lend themselves to use in localiza- 
БОП algorithms. 

im the prosecution of a submarine target, an air ASW 
crew will proceed through a succession of phases [Ref. 1]. 
The focus of this thesis will be the preflight, search and 
localization phases. The programs described in this thesis 
are designed to help the TACCO as well as Antisubmarine 
Warfare Operations Center (ASWOC) personnel make decisions 
regarding the selection of an initial search pattern. 
Additionally, a localization program is provided that could 
be used by the TACCO to expand a search pattern in order 


to convert from a single to multiple buoy contact. 
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The first two programs are designed to aid in pre- 
Flight planning. They transform the output from an acoustic 
model, such as ASRAPS III, FACT or PL-Ray, to a probability 
of detection plot and a graphical representation of the 
effectiveness of the buoy pattern selected. The graph gives 
the operator an indication of how well the pattern he has 
selected will work for his particular scenario. 

The remaining programs are used to aid the TACCO in 
the expansion of his pattern. The programs are designed 
to be used on an aircraft such as the P3-C (mod). Although 
Te is currently no provision for on-board processing of 
this nature on current aircraft, the upgraded computer 
systems that are in development should allow the addition 
of similar software. 

Іп general, the patterns employed will fall into three 
categories: barrier, distributed field, and convergence 
zone. The barrier pattern is designed primarily to gain 
direct path contact for a transiting submarine, while the 
distributed field and convergence zone pattern may be 
effective agains: either a transiting or on-station 
submarine. 

In barrier tactics a straight (or curved) line of 
buoys is employed at a position that is determined to be 
ahead of the transiting submarine. The pattern is 


monitored, and buoys may be replaced as they expire until 


ТЕГ 





Contact 15 gained or until it is determined that contact 
was missed [Ref. 2]. 
Advantages of barrier tactics include: 
1. The pattern is easy and quick to lay. 
2. Relative buoy positions are accurate. 


3. There is a possibility of initial multiple buoy 
contact. 


The disadvantages of barrier tactics include: 


1. The effective coverage areas of sonobuoys overlap, 
which reduces the total area searched. 


2. They generally allow only one pass of the target 
through the pattern. 


A distributed field pattern may be more desirable for 
an on-station submarine if its course and speed are not 
known and the pattern is based on an area searched. 

The general advantages given by distributed field 
patterns: 

i. The area of coverage is large. 
2. There is less overlap of effective coverage areas. 


3. There can be more than one chance of detection as a 
submarine transits the area. 


Disadvantages of distributed field tactics are: 
1. The field is ditficult to lay accurately. 
2. A long time is required to deploy the field. 


3. In a direct path environment, the pattern allows for 
less course variation than a comparable barrier pattern. 





Convergence zone tactics are also based on area Coverage. 
They may be either linear or planar arrays of buoys. Іп 
either case, they are based on gaining contact at a range 
distant from the sonobuoy. 

Advantages of convergence zone patterns: 

1. The area of coverage is large. 

A high probability of initial detection is obtained. 

The disadvantages of convergence zone patterns include: 

Meee Localization is difficult. 

2. The effectiveness of the patterns are dependent on 
the environment for existence of reliable convergence 
Zones. 

Additional scenario-dependent intormatior which must be 
evaluated includes the environment and geography of the 
area, and physical features such as straits cr landmasses 
Een may affect the prosecution of the target. Іп апу 
case, a knowledge of the environment is essential to a 
timely acoustic localization. 

The propagation loss curve is a primary planning tool 
for the TACCO. Combined with information about the target 
and sensors, it can be used to develop sonobuoy pattern 
spacings. There are some inherent problems in using the 
там propagation loss curve that arise because there is 
uncertainty in estimating the figure of merit (FOM). 

Actual detection ranges may differ from those forecast 


due to variations in any or all of the following factors: 


м 





1. The choice of acoustic model. 

Depth of the target. 

3. Actual frequency of target. 

4. Actual figure of merit. 

5. Beam pattern of the target. 

6. Other acoustic sources in the area of interest. 
;. Accuracy of environmental data. 

In a tactical situation, given a sound speed profile 
of the area and with a recent contact history, a search 
area can be reduced in size. This reduction of area is 
based on the concept that although the absolute magnitude 
of the figure of merit may not be known accurately, nor 
the actual levels of propagation loss, the acoustic models 
can give reasonably accurate range information. 

While it is obvious that the final objective is to 
track and attack submarines, a fact which normally requires 
multiple-buoy contact, a single-buoy contact may be converted 
to multiple-buoy contact by proper expansion around the 
buoy in contact. In the event that initial contact is gained 
on two buoys, much of the localization problem may already 
be soived. However, setting up a pattern for multiple-buoy 
initial contact requires many overlapping buoys which 
decreases the total area covered by the pattern as a whole. 
For a search pattern to be effective, multiple-buoy contact 


is not absolutely necessary, since the crew has the means 
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to deploy expansion buoys. Непсе, the measure of effective- 
ness used in this study will be the probability cf detection 
on a single buoy. 

The objective of this thesis is to provide a map of 
probabilities of detection associated with target position 
based on estimated target, environmental and sensor 
Barameters. 

The output of several of the programs are in the form 
of topographic maps which are based on a 40 by 40 rectangular 
grid. The programs allow for variable scales to be used. 
Mime sscale is defined as relative to the width of the map 
in nautical miles. The calculations required always use 
a 40 by 40 grid, allowing for finer resolution when the 


scale is reduced. 
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орки се LON MODELS 


Passive acoustic models use sound speed profiles to 
determine how sound energy radiating from an acoustic source 
is horizontally and vertically dis:ributed within the water 
column. The models also account for target frequency, 
target and receiver depths, bottom type and sea state. 
Additionally some models account for semi-coherent or in- 
coherent summing. 

A primary output of passive acoustic models is the 
propagation loss profile, which shows that the signal 


Intensity increases or decreases with range. 


РЕ ФФОКТЕ СИТТЕН DETECTION MODEL 

A basic detection model may be constructed using Figure 2.1. 
By calculating the FOM [Ref. 3], and drawing a horizontal 
line at that value, the ranges at which PL > FOM are those 
Where the propagation loss profile is above the FOM line. 
At these ranges, the signal excess is positive and probability 
of detection is unity (POD = i}. At other ranges, the 
signal excess is negative and the signal may not be discerned 
from the noise, hence, the probability of detection is zero. 
This type of model is called a "cookie cutter detection 
System". With it the effective area of a sensor has very 


well defined boundaries. 
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Bremer 22 2 shows the detection probability as a function 
of range for a cookie-cutter detector. This model is an 
idealized description of a real system. There are errors 
made in determining the FOM, which limits the accuracy of 
this model. The effects of inducing errors in the calcula- 
ЛОШ РОМ ате іп а large part dependent on the shape of 
the propagation loss curve. If the slope of the propagation 
loss curve is great near the intersection of the FOM line, 
an error of several dB in FOM will not greatly change the 
effective range of the sonobuoy. If however, the propagation 
loss curve is flat in the region of intersection, the 
effective range of the buoy may be markedly different than 


predicted. 


КОО RUSSIAN SIGNAL EXCESS DETECTION MODEL 

An improvement to the cookie cutter model is based on 
the uncertainty in predicted signal excess values. The 
Single look probability of detection is calculated from the 
predicted signal excess and sigma, the standard deviation 
of error in its estimation. The Gaussian Signal Excess 
Model [Ref. 4] assumes the error is normally distributed. 
The predicted signal excess is calculated from the sonar 
equation and a value for sigma is chosen, usually between 
mead 10 dB. Reference 5 gives values for typical systems. 

By estimating a figure of merit, a signal excess and 


a value for sigma, a probability of detection plot can be 
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created. This can be done by comparing the FOM and propagation 
loss curves at each range to obtain the average signal excess. 
The average sıgnal excess is then divided by the standard 
deviation to provide an input value x for computation of 
Ré Cumulative normal distribution function, ф(х), ѕіпсе іп 
the Gaussian signal excess model, the probability of 
detection (POD) = (x). 

The cumulative normal distribution function can be 


ExUnoxamated by the following expression from Reference 6. 


HEN) -5 (1€6€ X CcX'cX'«cx)" 


where 


C1=0.196854, 22=0.115194, C2=0.00034 and C4=0.019527. 


Once a plot such as that shown in Figure 2.3 is made, 
Ш сап be adapted for use in a number of ways. As a pre- 
flight tool, 1t will allow a ТАССО to visualize the effective- 
Mess Of a pattern of deployed sonobuoys. For example, a 
comparison between the plots for two buoy patterns could 
һе1р (0 determine which pattern is the more suitable, or 
it could indicate whether a pattern allowing small gaps 
between sonobuoys is preferable to one with no gaps but a 
Smaller area of coverage. 

It should be noted that estimates of FOM may vary 
during the course of a flight. It may be possible for a 


TACCO to improve the estimated FOM while in the search 
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phase, for example, if the ambient ncise meter readings update 
the value of the ambient noise level. The resulting change 
in the estimated FOM would also reduce value of sigma. 

For localization of a target to be effective, the value 
of sigma should initially be set as low as can be justified. 
Setting the value of sigma too high increases the area of 
ambiguity. While changing the value ot sigma alters the 
probability of detection at most ranges, ic does not affect 
the 50% probability of detection contour (MDR) on which 


buoy spacing may be based. 
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IIL. ORGANIZATION AND OUTPUT OF PROGRAMS 


The programs developed for this thesis are designed to 
give graphic outputs, which may be quickly but accurately 
interpreted зу a trained operator. Graphics are desirable 
because they can convey much information in a short time. 

The organization of this analysis makes possible the use of 
compatible data files to store data that will be calculated 
by one program and transferred to another. The use of text 
files ailows this storage with a minimum of errors and 
operator input time. The use of disk storage allows a 
proragation loss curve to be manually entered once. It is 
then retrievable for any analysis including a variety of 

maps and graphs comparing different propagation loss profiles. 

Figure 3.1 shows the operational layout of the programs 
and data files. The inputs to the acoustic model are estimated 
as in the preparation for a normal preflight briefing. The 
accustic model may be any that provides a propagation loss 
Enortrle. 

There are two text file types used and they may be 
i ur ic Бу the suffix "-PLF or -HIRES". The suffix is 
a code identifying the type of data stored in the file. 

It should be noted that each text file is merely a form 
of data storage and retrieval and does not operate mathematically 


om the data. 
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A.  PLF-MAKER PROGRAM 


D rst program 1s PEF-MAKER. It prompts the user for 


a filename and for propagation loss data at one nautical 
mile increments for ranges of 0 to 50 nm. The operator 
tne option of correcting the profile prior to final 
storage. Їп this analysis the riles are coded as follows: 
ATL 050 5 р 5 
NNN fff Td Hd s -PLF 
NNN - Name of ocean area 
ATL - Atlantic 
BDA - Bermuda 
ШОЕ = Indian Ocean 
MED - Mediterranean 
PAC - Pacific 
IS Frequency 
0509 7505112 
Td - Target Depth, S = Shallow D = Deep 


Hd - Receiver Depth, S = Shallow Р = Deep 


s - Summing, s = semi-coherent 
1 = incoherent 
PLF Propagation Loss fiie 


HIRES - High Resolution Plot Datafile 


B. PLOT I PROGRAM: 
The Plot I program prompts the operator for inputs for 

two cases. The inputs are propagation loss "-PLF" files, 

figures of merit and sigmas. The program provides a 


graphical plot of the propagation loss curve for each 
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meee This is a graphical plot of the data provided to the 
PLF-maker program. 
ги сссопа разїт of plots is created by calculating: 
P (SE > 0) 
for each range increment, which is based on the estimated 
FOM, sigma and the above propagation loss profile. 

A third plot is a difference of probability of detection 
plot. This may be seen as an improvement plot. It is a 
point by point comparison of two percentage of detection 
plots. Positive values indicate that the second case has 
a higher probability of detection at a given range and a 
negative value indicates the first plot is more effective 
at that range. This plot is useful because it shows the 
TACCO the areas where improvement or degradation may be 
expected, if the parameters calculated in the FOM may be 
updated. 

Figure 3.2 shows the output from the PLOT I program. 
In this case it demonstrates the effect of increasing the 
ute Of merit from 75 to 78 dB. This may occur due to 
the TACCO receiving updated environmental information, 
changing processing modes, or both. If updated ambient 
noise values are obtained, the FOM and value of sigma may 
both change. 

This program demonstrates the effects of changing any 


or all of the parameters used to calculate the probability 
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Figure 3.2. Output from PLOT I Program 


First Propagation Loss Profile 
Second Propagation Loss Profile 
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of detection plot. The program allows the operator to compare 
two propagation loss files, for example, different models 

or different frequencies using the same model. Two figures 

of merit may be compared, or two values of sigma may be 

used. The results of each of these changes may be computed 
 Еапасеіу оу they may be taken as a group. This program 

is designed to answer "what if...''-type questions for the 
operator. It can heip him determine those areas in which 

an error in estimating figure of merit has the greatest 
effect. 

The program may also demonstrate the effects of varying 
the value of sigma. If the value of sigma is increased, 
the probability of detection function flattens out and, 
in an extreme case, approaches 50% for all ranges, giving 
an ambiguous position. Target contact provides very little 
meealization information. 

While the value of sigma is generally taken to be 5 to 
10 dB, this program demonstrates the effect of increasing 
or decreasing sigma. As sigma is decreased, the percentage 
detection approaches a step function or "cookie cutter" 
detection system and contact on a single sensor produces 
a relatively small area of possible locations for the 
ameet. Іп а "real world" situation this could occur if 
all attributes of the target and environment were well 
known and if temporal and spatial variations were non- 


existent. Most often, the cookie cutter model is 
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Pea, пос because of its accuracy, but because of its 
Simplicity. 

A practical use of PLOT I is to compare deep and 
shallow hydrophones with a given target depth and sound 
Speed profile. This may be used to relate to the operator 
the advantage of using one or another hydrophone depth 
tins. 

Two values of FOM may be input for the same curve 
showing the importance of correctly calculating this 
Parameter. This type of analysis can show the operator 
how much increase or decrease in range he can expect due 
to changing integration times on the acoustic processor 


[Ref. 7]. 


wee MAP SERIES 

Mme output of the MAP series of programs are data 
mbes Containing probability of detection for an area of 
pic етаїһ'<$ surface. A square shaped area is divided into 
a 40 by 40 grid, thus giving 1600 cells. The size of the 
S 1S determined by the operator input of scale, which 
is the length of a side. Generally the scale used will 
be 50 nm, although the programs allow the operator to 
input any value. If the scale is increased greatly, some 
resolution may be lost, but in a localization scenario, 


scales may be successively reduced for better resolution. 





The programs use buoy positions in x,y coordinates 
Ben values in nautical miles. The center of the plot is 
tor all scales. Shown in Figure 3.3. are x and y 


values for the 50 nm scale. 


D. THE MAP 0 PROGRAM 

The MAP 0 Program provides a graphical anaylsis of a 
buoy pattern. The output plot may be used to show the 
effective extent of a sonobuoy pattern. it also shows 
how the detection probability varies throughout the area. 
The highest single buoy percentage detection is plotted, 
hence the program assumes complete independence among the 
sonobuoys. 

The output is a 'HIRES" datafile which is compatible 
with the CONTOUR II and BLOCK II programs as modified for 
this project [Ref. 8]. When used with CONTOUR II, the 
program MAP 0 gives an output of target coverage, which is 
designed to be useful in mission planning. The pattern 
drawn indicates areas of high and low probability of contact. 
The density of the contours indicates the accuracy of 
phe prediction, with closer packing indicating greater 
accuracy. This may be used to show probability of detection 
history if the target track is plotted as an overlay on 
the map. 

The isopleths, as shown in Figure 3.4, are lines of 


equal probability; generally lines will be plotted at 
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ке 5, CONTOUR II Plot of a Single Sonobuoy 
Detection Probability 


20) 





equal intervals. In instances where the slope of the 
percentage detection plot is very steep, intermediate 
pues may be eliminated for clarity. 

In the example shown, the 50$ contour is a continuous 
line indicating that a target may not penetrate the pattern 
without being exposed to. 0.5 probability of contact. This 
program may be used to give a plot of the target coverage 
of a complex pattern of sonobuoys. Chapter 5 includes 
analysis of patterns developed by NADC to test a localization 
algorithm [Ref. 9]. Various propagation loss profiles are 
compared for the different patterns. 

Figures 5.4 and 3.5 are contour and block representations 
of probability of detection for a single sonobuoy. The 
horizontal scale of each of the displays is 50 nm. Figure 5.4 
shows two bands within which the probability of detection 
is greater than 25$ located at 12 and 20 nm. The block 
diagram, Figure 3.5, shows that the rise in probability 
is slight at those ranges. The probability of detection 
plot used for this example is from Figure 2.2 

An additional use of the MAP 0 plot is in planning 
Simulations. If the known track of the target is plotted 
directly over the contour chart, the probability of 
detection may be determined at any time during the run. 
Single sonobuoy plots may also be made and a determination 
of the probability of contact on the lower priority sono- 


buoys may be determined. 








amure 3S BIOCK II Plot of a Single Sonobuoy 
De heelroner Toba Bl tty 





E. MAP I PROGRAM 

LENS the development Of a Localization algorithm 
Eus the propagation loss profile and contact status for 
determination of position. Any number of sonobuoys may be 
used and any number of those may be holding contact.  DIFAR 
information is not used by the algorithm. 

The algorithm calculates the probability that the buoy 
status is as observed assuming the target is in the center 
of each given cell. Positive contact on a buoy causes a 
peak to be formed near the buoy while negative contact on 
another buoy causes a depressed area to form in the vicinity 
of that buoy. 

It should be noted that in general the probability of 
detection at long ranges is near 0, except in convergence 
zone environments. It may be possible for the operator to 
delete buoys that clearly have low probability of contact 
to decrease the run time of the algorithm. These buoys 
add little information about the possible location of the 
target. 


The probability calculated is equal to the following 


product: 
n m 
J{ POD (buoy in contact) ~ 7{ 1-PoD (buoy Ji Contact) 
1 ШЕ! 


where n = number of buoys in contact and m = number of buoys 


їп pattern. 
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The display will indicate the areas to concentrate on 
as potential sites for additional sonobuoys. The map thet 
is output shows peaks and valleys. The peaks indicate 
areas Of high relative probability of contact. Thea Acco 
may elect to deploy a sonobuoy on one of the peaks in which 
case, two outcomes are possible for the new (n+l)st buoy. 
If contact is gained on this new buoy, then POD(n+l) becomes 
(Гог іп the above expression. Since this factor is 
close to 1 for short ranges and close to 0 for distant 
ranges (assuming no convergence zones), the fact that 
contact is gained has the effect of decreasing the height 
Of any distant peaks that may exist. If no contact is 
gained, 1-POD(n+l) becomes a factor. This factor has the 
r to effect. The local peak is decreased and distant 
peaks are not affected. In the existing program it is 
desirable to eliminate buoys which do not hold contact 
and which are clearly out of range of the target. 

The resultant HIRES file, when plotted by CONTOUR II 
or BLOCK II, denotes the areas most likely to contain 
the target for an existing pattern of contact status to 
arise. Peaks will arise as closed contours. The areas 
indicated by this map should then be searched by expansion 
buoys to complete the localization of the target. The 
Becuracy Of such a localization algorithm is dependent on 


how well the acoustic model duplicates the environment. 


54 





Single buoy DIFAR contact may be used to aid the TACCO 
by indicating which of several peaks are most likely to 
contain the target. However, in the event of unstable or 
mecurate bearings, the DIFAR information could hide an 
otherwise existent peak. In the case where the environment 
is improperly modeled, this method may not agree with the 
DIFAR information and confuse a merged routine. In the 
event that DIFAR contact is held on multiple sonobuoys, 

a Kalman based localization algorithm such as ASW MUSCLE 


has been shown to be effective [Refs. 9 and 10]. 
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IV. ANALYSES OF VARIOUS ENVIRONMENTS USING PLOT I 


This section describes the origin of the data base used 
dan terprets the output of the PLOT I program. Various 
acoustic models were used to demonstrate the flexibility 
of the program. The comparisons that follow are made based 
on the propagation loss curves using the data in Example I 
combined with sound speed profiles characteristic of the 
areas shown. The first comparison is a FACT model propagation 
loss curve. PLOT I is used to demonstrate the effects that 


changing the various parameters have on detection probability. 


EXAMPLE I. Improvement of Figure of Merit 


Filename PASOS 
Ocean Area к Pacific 
Model : ASRAPS III 
Latitude 35 50 N 
Longitude 240125: N 
sea State : 2 
Bottom Province 2 3 
Sonic Layer Depth : ШОКЕ Е 
Target Depth E ВОО ft 
Hydrophone Depth : 90 ft 





Figures 4.1 and 4.2 are conventional propagation loss 
curves. The horizontal line represents the caiculated 
ieee Of merit. Currently this is the only type of plot 
avallable to the TACCO during the preflight phase. The first 
plot shows an MDR of 7 nm and a possibility of convergence 
zone detection between 35 and 38 nm. The propagation loss 
curve rises at that range but it is impossible to determine 
the likelihood of gaining contact beyond stating that the 
probability of detection is less than 0.5 

Figure 4.2 shows that by increasing the FOM the MDR is 
increased to 10 nm and an improvement in probability of 
detection (nearly 0.5), but since the POD was not able to 
be calculated in the previous case, the amount of improvement 
Still cannot be determined. 

Together these plots can be used to show the improvement 
that may be expected by increasing the FOM from 75 to 78 
dB. Determination of the MDR may be made by extending a 
horizontal line at the 50% probability of detection (see 
Figures 4.3 and 4.4). This method will also determine the 
POD in convergence zones. For the direct path portion of 
the curve, the greatest increase occurs near the MDR. 
ае 4.3 shows а 0.20 to 0.25 probability of detection in 
the convergence zone. By increasing the FOM (Figure 4.4) 
the POD increases to nearly 50%, a significant increase. 

КОО e 45 15 a "Difference Plot". It is a comparison 


of two percentage detection curves (in this case Figure 4.3 
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and 4.4). It shows where the two curves are in agreement and 
where they disagree. This may be interpreted as an "improve- 
ment" plot wherein the second data set is compared to the 
rst. In cases where the plot falls below the zero line, 
the first plot is more likely to gain contact than the 


second. 


КОРЕ ТТ. Improvement of Figure of Merit 


Filename : 10005510) 015 
Осеап Атеа Indian Ocean 
Acoustic Model : РІ,-Кау 
Latitude : 03 00 М 

re; State : 3 

ton Province ; 3 

Sonic Layer Depth : JEONG) ET 
Target Depth 500 ft 
Hydrophone Depth =: 1000 ft 


Figures 4.6 and 4.7 show a typical Indian Ocean sound 
speed profile using the data from Example II. The FOM's 
analysed are 75 and 80 dB. Based on the propagation loss 
profile, the operator can see that there is a possibility 
of bottom interaction propagation in both cases. There 15 
а greater probability in the second case but the magnitude 


of the probability may not be determined. 
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The next step, Figures 4.8 and 4.9, show that the second 
case (80 dB) has a greater possibility of detection for all 
ranges. The areas of greatest improvement occur in the 
regions of peaks in the propagation loss ae The 
probability of contact is improved from less than 25 percent 
to about 50 percent. 

Figure 4.10 shows that a large improvement over most 
ranges occurs when the FOM is increased from 75 to 80 dB. 

By comparing Figures 4.8 and 4.9 the MDR does not change 
appreciably, it stays in the 2 to 53 nm range. However, in 
Ene latter case the plot rises and dips greatly, indicating 
a greater possibility of intermittent or weak contact if 
the target transits at ranges from 6 to 22 miles from the 
buoy. 

Figure 4.11 is the same propagation loss profile shown 
ши Езситез 4.6 and 4.7, but with a FOM of 85 dB. This 5 dB 
change increases probability of detection through all ranges 
Em Uown by Figure 4.12. There is now a 75 percent 
probability of detection in the bottom bounce regions, and 
there is also an improved probability of detection at all 
ranges. Figure 4.13 shows that in this case the improvement 
from 80 to 85 is of more benefit than that from 75 to 80 
for this region. This would give the mission planner an 
Enuscatyon of the advisability of using a search pattern 
based on bottom reflected paths. If he is relatively 


sure his FOM will be at least as good as 85 dB, he may 


42 





50 


60 


70 


ав 
90 


90 


100 


50 


60 


70 


ав 
90 


90 


100 


Эн 


10 20 30 40 20 
RANGE 
igure 4.0 Propagation Loss Profile 
106050 5 FOM: 75 Sigma: 5 
| 
10 20 30 40 50 


КАХСЕ 


Piemmes4.7 Propagation Loss Profile 


ТОСО5Орр5 FOM: 80 Siema На 


45 








75 


POD 


.50 


225 


Figure 4.8. 


iu ¥ 20 


\ 


| | 
10 20 


Figure 4.9. 


30 
RANGE 


م س 
30 
RANGE‏ 


ЕОМ: 


40 


75 


40 


Probability of Detection 
IOC050DDS 


50 


Sigma: 5 


20 


Probability of Detection 


IOCOSODDS 


44 


FOM: 


80 


Sigma: 


5 





200 


025 


-0.5 


MA m а 


10 20 30 40 
RANGE 


те 1.10. Difference Plot 


Case 1 IOC300DDS FOM: 25 
Case 2  IOC300DDS FOM: 80 


45 


Sigma: 
штоша: 


20 





50 


60 
70 
dB 
80 
{N\A : 
y N f AAN 
90 
100 | 
“0 10 20 30 40 50 
RANGE | 
шэн. osse Profe 


IOCO5S0DDS РОМ к КЫз Улы ста: 5 


‚PoD 





225 


0 : 
0 10 20 30 40 50 
RANGE 


Figure 4.12. Probability of Detection 
IOCO50DDS FOM: 85 Sigma 5 


46 





2-0 


05 
E VV rn 
-0.5 ` 
0 
0 10 20 30 40 50 


RANGE 


Бирите 4:15, “Difference Plot 
Casei: TOC300DDS FOM: 80 Sigma: 5 
Case 2 LOCS00RDS БОГ, Sigma: 5 





eT to use a larger pattern spacing to increase the total 
area searched. If he is not confident of the FOM then he 
Should use a direct path pattern although the area coverage 


is less. 


EXAMPLE III. Variation of Sigma 

This section uses the same propagation loss profile as 
Figure 4.6. The two values of sigma used are 10 and 5 dB, 
the results being shown in Figures 4.14 and 4.15, respectively. 
It should be noted that when SE = 0 there is no effect on 
the probability of detection. The effects noted in decreas- 
ШІБ the standard deviation is to increase probability of 
detection in areas of positive SE (0 to 3 nm in this example) 
КІ to decrease the probability of detection in areas of 
negative SE (all ranges greater than 3 nm). When SE has a 
large absolute value (greater than 2 times the larger sigma) 
There 1s little or no affect in decreasing sigma. 

In Figures 4.14, 4.15 and 4.16, the same propagation 
loss profile was used, with only the standard deviation of 
the FOM changing between the two ceses. The difference plot 
shows an apparent degradation of detection of nearly 20 percent 
from 5 to 30 nm and a smaller effect at longer ranges. 
Since decreasing the standard deviation is the same as using 
a more accurate model, lower probabilities should not be 
interpreted so much as a "degradation" in the detection 
capability of the sensor system, but as an increase in the 


ability of the system to localize, once contact is gained. 
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The difference plot is useful because it identifies the 
areas where knowledge of the environment has the greatest 
effect on localization. In this case decreasing sigma from 
Шо 5 decreases the POD about 20 percent, except at ranges 
from 25 to 50 nm. This is due to a depression in the 


ition loss profile at these ranges. 


EXAMPLE IV. Variation of Hydrophone Depth 

This case is based on two different propagation loss 
profiles for the same location. The same model and input 
mrOormation are used for both runs with the exception of 


hydrophone depth. 


Filename BDAOSOSDS BDAOSODDS 
Ocean Area Bermuda Bermuda 
Model DE Ray PL-Ray 
Patıtude 36.00N 36.00N 
Bottom Province 4 4 

Шлгрес Depth 400 ft 400 ft 
Hydrophone Depth 90 ft 1000 ft 


Figures 4.17 and 4.18 show the propagation loss profiles 
for both cases. The deep hydrophone case has a longer 
MDR but it drops more rapidly and fluctuates greatly as 
range increases. The shallow hydrophone example exhibits 
a smoother behavior throughout the range. If one inspects 
fie propagation loss curves alone, it is difficult to tell 


if either depth setting is preferable over all ranges. 


Ed 





Using a FOM of 75 and a sigma of 5 results in the 
Морава: тту о detection pröfiles in 4.19 and 4.20. Both 
are very similar with the exception of a longer MDR and 
broader CZ for the deep-deep combination. 

Figure 4.21 confirms the greatest change difference 
between the two plots occurs at the short range where the 
difference is large. At the CZ range, although the deep 
hydrophone has a broader annulus, the probability of 
detection is higher on the shallow hydrophone for a short 
distance. By comparing the percentage of detection plots, 
the shallow hydrophone shows a maximum of 0.3 POD at the 


MAN le the deep shows only a 0.25 POD. 
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V. COMPARISON OF PATTERN EFFECTIVENESS 


This chapter demonstrates the utility of the MAP 0 
Program in Simulation and in flight planning. The patterns 
and tracks used were taken from Reference 8 and are typical 
of those used in initial search. Various propagation loss 
profiles are used to show hew the effectiveness of the 
patterns are affected by the environment. 

Reference 9 showed that certain patterns exhibited 
variable success in detection and localization. This may 
be attributed to the fact that the contact effectiveness 
of a pattern is quite dependent on the environment that is 
used. In simulation planning the MAP 0 Program ailows the 
initial search pattern and environment to be evaluated prior 
to the run, by superimposing target tracks directly over 
the MAP 0 output. 

The reader is instructed to refer to Figures 4.1, 4.3, 
4.6 and 4.8 for the propagation loss curves and percentage 
detection plots used to develop the MAP 0 outputs. The 
PAC300DSI (Pacific Ocean) environment will be used to 
demonstrate MAP 0 and outputs for appropriate patterns. 
The same patterns using the IOC300DDS (Indian Ocean) case 
will be used to demonstrate how the pattern may become 


inadequate if an alternative environment is used. 
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Figure 5.1 shows three representative target tracks 
relative to eight sonobuoy positions. This plot was made 
using the PAC300DSI propagation loss profile. For simulation 
planning, a buoy pattern may be input with the required 
FOM, Sigma and Propagation loss Curve. 

The first environment (Figure 5.1) makes an efficient 
pattern, suitable for both detection and localization. 
Target #3 enters the pattern from the East and the probability 
or detection rises until CPA on the first sonobuoy. The 
50% contour is an unbroken line and the 70% contour is 
large, although broken, the target may not move within the 
pattern without being exposed to at least 50% POD. 

In the second envirunment (Figure 5.2) it may be noted 
that target #3 has a very high probability of detection on 
the first buoy along his path but that the target never 
quite reaches the 50% contour on the second buoy before 
firming. This probability of gaining a detection on the 
third buoy is even lower. The pattern is clearly not 
adequate for this environment. The gaps are large and it 
is likely that a target could enter the pattern and be 
missed entirely. 

Using MAP 0 and visualizing possible paths through the 
pattern, the TACCO may predict the effectiveness of the 
sonobuoy pattern and identify any weaknesses. By changing 


buc spacing, the TACCO can change the effectiveness 
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of the patterns. The crew is limited in the number of buoys 
that can be monitored effectively and decreasing the spacing 
decreases the total area that can be searched. Using this 
map and estimating the Ши гс position or probability 
area, the TACCO has a graphical display cf the effectiveness 
of his pattern and any gaps in coverage that may exist. 

A typical barrier pattern is shown in Figure 5.3 The 
plot shows that in this case there are "gaps" in coverage 
between the buoys where the probabiliby of detection falls 
below 50%. This may be remedied by decreasing the spacing 
between buoys, however, that will also decrease the width 
ШШ Пе search area. In flight planning the apriori distri- 
bution of the target may be used to determine if tle pattern 
Should be compressed. 

TT the apriori distribution is large compared to the 
buoy pattern, the gaps that exist may vary in size and 
extent. The example in Figure 5.3 has gaps that are 
relatively narrow, only one tenth of the bucy spacing, and 
relatively unimportant, the probability of detection is 
40% at its worst. The example in Figure 5.4 shows the 
Same spacing, but in this environment, the pattern is much 
less effective, there are large gaps and a target penetrating 


the barrier may see a probability of detection as low as 10%. 
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VI. DEVELOPMENT OF MAP I, A LOCALIZATION AID 


iE output of the MAP I Program is a topographic 
representation of the probability of detection calculated 
over the surface. The program displays a surface based on 
either pcsitive (C) and negative (N/C) contact or positive 
information alone. The input is the propagation loss profile, 
FOM, sigma, buoy pattern and contact status. 

For each of the 1600 cells, the program computes a 
probability and assigns the value to the respective cell. 
A representative 3 buoy pattern would allow the surfaces 


based on the possible contact status combinations as shown 


below: 
Buoy Status Parameter Calculated 
#1 #2 #3 
E N/C N/C PML) (1 POD(2))* (1 PORK 3) ) 
G C М/С БО (1: 15200(21"(1-РОМ В 77) 
S E C EDIL) SPODE2)*POD(3) 


The single look probability of detection is calculated 
as if the target were located in each cell and the resultant 
Шишрет assigned to the cell. The probability surface is 
plotted. Peaks on this surface are indicative of areas of 
high relative probability of detection. 

It was noted that on some propagation loss profiles 


values may change by 20 dB over a short range. To avoid 
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missing a peak the probability of detection as a function 
of range is replaced by an envelope function which takes the 
greatest of 3 POD computations for each range. 


That is, for each rane, r: 


Max(POD(r-1 nm) ,POD(r).POD(4+1 nm)) 


where r is in nautical miles 


ENcalculated for the probability of detection function for 
ws in contact. A similar function is calculated for buoys 


Hol ІП contact: 


Max(1-POD(r-1 nm),1-POD(r),1-POD(r+1 nm)) 


The result of using the envelope functions are larger 
probability areas, but they decrease the problem of missing 
detections due to a greatly varying environment or a rapidly 
fluctuating propagation loss curve. Although it increases 
the area that remains under suspicion, it reduces the effects 
of underestimating the propagation 1oss curve which could 
lead to missing an area of high detection probability. This 
method also decreases the problem of inaccuracies in range 
Si the propagation loss curve. 

The program has the ability to use positive and negative 
information, that is, buoys holding contact will increase 
the value of the probability in the local vicinity and 
decrease it at longer ranges. A buoy not holding contact 


will decrease this parameter primarily at close ranges. Тһе 
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effect at long ranges will be much less. If two buoys are 
in close proximity with the first holding contact, the surface 
will peak near the first but will be skewed away from the 
second. Figure 6.1 shows such a case. This method seems 
to be most effective when only one or two buoys are holding 
contact. 

A second 11:27 CoO accEeDE Only positive data. 
This is less accurate but run time is decreased. It may only 
be effective when two or more buoys are holding contact. In 
such a case, the absolute values of the entire surface may 
decrease drastically because most of the probabilities are 
much less than 1, and when multiplied decrease greatly. 
Peaks still occur where high probability areas intersect. 
Figure 6.2 shows the same pattern as above, with both buoys 
holding contact. The peak has now shifted from outside of 
the two buoys to the center of the two buoys. This information 
may be used by the TACCO to plan expansion buoys. If DIFAR 
information is available, it should decrease the size of the 
probability area even further. 

The following development shows the operation of the 
MAP I program as a target passes through a buoy field. 
The target and pattern (Figure 6.3) are taken from Reference 
10, гип 17, It should be noted that the output of the 
program is predictable from the contact status alone and is 
not dependent on the observed signal excess. If, for example, 


the FOM were less than predicted, the contact status would 








Figure 6.1. MAP 1 Plot, Two Sonobuoy Pattern. 
buo! Holbldrnescontact. 





кете 0220 MAP I: Two Buoy Contact 
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М ре at different times, but as long as the buoys gained 
contact, the MAP І output would be the same. 

The track used here is the same as track 3 in Figure 
2.53. AS the target enters the buoy field, initial contact 
is gained on buoy #3. Figure 6.4 shows the MAP 1 output 
when contact is initially gained. The peak is skewed to 
the east, away from the rest of the pattern, since contact 
is held only on one buoy. 

As the target proceeds along its westward track, contact 
1s gained on buoy #2. Figure 6.5 shows the output of MAP I 
with this new information. The peak is between the two 
Semeact buoys but is slightly north of the target's track. 
In this case the peak is close enough to the target to give 
adequate localization. 

After the target turns, contact is gained on buoy 5. 
This causes the peak to shift to the south and is positioned 
roughly in the center of the three contact buoys (Figure 6.6). 
Again the peak is near the track of the target. 

As the target continues, contact may be held briefly 
on four buoys. In that case the output would look like 
Peeure 6.7. In the event contact was lost on any one of 
the four, the peak would be moved away from the individual 
sonobuoy. 

The MAP I program could be used in localization to help 
the TACCO determine where to lay the next sonobuoy when 


Banane contact is gained. Given the current buoy contact 
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status, he may add an additional buoy at or near the peak. 
Mims DUCy has two possibilities; it either will gain contact 
Or 1t will not. In the first case Sie peak will be slowed 
moma position between the buoys in contact. If contact is 
not gained, the buoy peak will be depressed and other lesser 
peaks will remain. The TACCO may then deploy sonobuoys in 
the vicinity of these until multiple buoy contact is 
achieved. The advantage of this program is that it gives 

the ability to proceed in a logical order to investigate 


areas of probability. 
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Figure 6.3. Sonobuoy Pattern and Track ot Target 





Паше 0.4 Contact Gain on Buoy #3 
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Figure 6.5. Contact Held on Buoys #2 and #3 


+ a; Me 
* * * 


шог 16 Contact Held on Buoys #2, #3 and #5 








Бисште 6.77 Соптасс Held on Buoys #2, #5, #5 and #8 
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VII. CONCLUSIGNS 


Various products can be created by using the signal 
excess model. The plots produced by the programs based on 
this model utilize microcomputer graphics to provide infor- 
mation which may be used to optimize sonobuoy employment. 
The inputs required are a propagation loss profile, a 
figure of merit and a standard deviation of the signal excess. 
The plots shown can indicate the effectiveness of a sonobuoy 
pattern. 
The localization program has several pctential advantages. 
These include the following: 
1. It is independent of DIFAR processing. 
2. It leads the aircrew to areas ot high probability. 


3. It may be used in direct path, bottom bounce and 
convergence zone environments. 


4. An inflight version could use bathythermographic data 
coupled with an acoustic model to produce more accurate 
results. 


5. It is flexible in that it may be used in any environment. 


There are also a number of limitations that are evident: 


1. It requires long run times. Machine level programming 
шин ер ап бретастбпа1 ртортат рїүе а геа1 time output. 


2. It requires much interaction to input data. Direct 
computer interfacing that currently exists could be 
used in an operational program. 


3. As with any localization program, false alarms and 


missed detections could have an adverse effect on the 
program. 
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There are also some considerations that could be 
ENcorporated to improve the efficiency of the localization 
program. Multiple propagation loss profiles would allow 
the correlation of different frequencies. A directional 
ambient noise model such as the DANES model would add some 
capability. However, DANES would increase the required run 
times, since it would need a much larger matrix to describe 


the directional percentage detection profile. 
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APPENDIX A 


The following program listings will provide the plots 
shown in this thesis. The programs are written in Applesoft 
BASIC and require an Apple II computer with 48K and a single 
Шек drive. For the BLOCK II and CONTOUR II programs only 
the changes required to convert them from the BLOCK and 


CONTOUR from Reference are listed. 


2 





та 
29 
за 
48 
5й 


79 
88 
20 


PLF-MAKER PROGRAM LISTING 


DIM PLi SO) 


D$ = CHRE (4) 


INPUT "MAME UF PRÜFLOSS CURVE FILE :"iN¢ 
FOR I - 9 TO 58 
PRINT “INPUT PROPLOSS AT ";l;" MM: "i: INPUT ""; 
PL X I) 
NEXT I 
PRINT ОФ"ОРЕМ "МВ + "-FLF" 
PRINT DS"WRITE "¡N3 + "-PLF" 

FOR I = 8 TO 58 

PRINT PLI? 

МЕХТ 1 

PRINT D$"CLOSE “srt + "-PLF" 

GOSUB 300 

INPUT "ANY CORRECTIONS (YM) 2":69: IF é$ = "N" 
ТНЕМ 200 

INPUT "RANGE : "iR 

INPUT "PROP LOSS VALUE : ";PL&R) 

INPUT "ANY MORE CORRECTIONS (YM) ? :";A$: IF A 
$ = "N" THEN 88 

ТЕ A$ = "Y" THEN 158 

END 

PRINT D$"OPEN ";N$ + "-PLF" 

PRINT D£"READ "Nt * "-PLF" 

FÜR I = 8 TO 58 

INPUT РЪСТ) 

NEXT I 

SPEED= 100 

FOR I = 8 TO 4? STEP 2 

PRINT I,PLCD ,PLCOI + 1) 

NEXT 1 

PRINT S8 ,PL{ 58) 

PRINT D$"CLOSE "N$ + "-PLF" 

SPEED= 255 

RETURN 
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18 R 
20 D 
110 


120 
140 


158 
158 


120 
200 
300 
310 
320 
334 
340 
330 
354 
379 
сай 


518 
520 
534 
550 


560 


665 
678 
гав 
214 
229 
4020 
4425 
4027 
4028 
4030 
49490 


44356 
4668 


4429 
4080 


PLOT I PROGRAM LISTING 


2 = 59 
Ф = CHRE (4) 

Ci = .198854:С2 = ,115194:С3 = .000344:С4 = .815 
227 | 

FORTI - 1 TO 2 

иг 0 МЕМЕЗБЕ ЕТСЕ Hosts’ 1:32 INPUT ""; 
MAS CID | 

БЕТЕ ЕЕЕ О sie INPUT “> 3" sFMCI> 

PRINT "INPUT STANDARD DEVIATION C":I;:: INPUT "2 


".8661) 
МЕХТ 1 
DIM PLC2,RZ2»? ,GC2, R2 ,DFCR2D 
FOR J = 1 TO 2 | 
PRINT D$"OPEN ";NAà$(J) * "-PLF" 
PRINT D$"REAÓD ":Né$OD + "-PLF" 
FOR I = 0 TO R2 
INPUT PL(J,D 
NEXT 1 
PRINT D$"CLOSE ";N@SCJ) + "-PLF" 
NEXT J 
PEE сс = j 7 <2 4 <1 + C1 Xx X + C2 x 
ТО тг s + C4 X CX * 4D) ^ 4 
FOR I = 8 TO 56 


FOR J пощите 
Bo FMCJ) = PLJ.IO 

Ja eee HM ИМЕН Ве. Те INT (188 X FN P 
SEM TSG SE JD: SOT 665 
Sout = ПМ C190 — (188 X ЕМ РФ - БЕС) / 55% 
u5 229 

NEXT J 

NEXT 


FOR I = @ TO R2 
РЕСІ? = 041,1) - 9‹2,1) 


МЕХТ 1 
FOR J = 1 TO 2 
HGR 
HCOLOR= 3 

PN = 1 


HPLOT 8,8 TO 8,1568 TO 279,150 TO 279,98 TO 8,8 

HPLOT 0,30 TO 5,30: HPLOT 275,38 TO 279,38: HPLOT 
9,60 TO 5,58: HPLOT 275,48 TO 279,60 

HPLOT 8,98 TO 5,96: HPLOT 275,98 TO 279,98: HPLOT 
8,120 TO 5,128: HPLOT 275,128 TO 279,129 

FOR I - 5 TO R2 STEP 5 

Оо 2 2 270 282,144 TO I x 279 / R2,149 

NEXT 
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4085 


4087 
4090 
4100 
4105 


4110 
4120 


4125 
4130 
4 1208 
4156 
4170 
4175 
4 189 
4215 
4229 
4225 
4230 
4240 
4250 
4259 


4261 
4252 
4253 
4270 
4280 
4270 
9300 
4345 
4306 
4310 
4315 
4429 
4436 
4440 
4459 
4454 


HPLOT B,(3 X FM(J) - 1580) TO 279,(3 % ЕМОр - 
150) 
IF PL¢J,@) < 59 THEN HPLOT 9,0: GOTO 415# 
HPLOT 0,3 X (PL(J,0) - 50) 
FOR I = 1 TO R2 
IF PL(J,1> < 58 ТНЕМ HPLOT I X 279 / R2,8: GOTO 
4130 
IF PL¢J,1)» > 106 THEN GOTO 4125 
HPLOT TO I X (279 / RD,3 X <PL{J,D - S8»: GOTO 
4 130 
HPLOT TO I X 279 / R2,150: GOTO 413A 
NEXT I 
PRINT 
PRINT "PROPLOSS CURVE "¡NAs<J) 
PRINT "FOM : ";FMXJD;" SIGMA : ":5б62) 
INPUT Hš: IF Hš = "P" THEN  GOSUB 59909 
NEXT J 
FOR J = | ТО 2 
HGR 
PN = 2 
HPLOT 9,8 TO 9,158 TO 279,154 TO 277,9 TO 4,4 
HPLOT 8,75 TO 19,75: HPLOT 269,75 TO 279,75 
HPLOT 8,27 TO 5,37: HPLOT 274,37 TU 277,37 
ELT 0,111 TO 5,111: HPLOT 274,111 TO 279,111 
FOR I = 5 TO R2 STEP 5 
HPLOT I X 279 / R2,144 TO I Х 279 / 82,149 
NEXT 
HPLOT 8,158 - 1.5 х 957,9) 
FOR I = 1 ТО Р2 
HPLOT TO I X (279 и Е2) „158 - 1.5 x (000,13 
NEXT I 
PRINT "PERCENTAGE DETECTION PLOT " ; мањ ‹ Ју 
PRINT "FOM : "sFMCJ> :" SIGMA : 5600 
INPUT HS: IF H$ = "P" THEM GOSUB 5000 
NEXT J 
HGR 
HPLOT 8,8 TO 8,150 TO 279,158 TO 27?,8 TO 2,8 
HPLOT 9,75 TO 279,75 
HPLOT 8,27 TO 5,37: HPLOT 274,37 TO 279,37 
HPLOT 9,111 TO 5,111: HPLOT 274,111 TO 279,111 
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4441 
4462 
4463 
4476 
4480 
4476 
4588 
4503 
4505 
4504 
4507 
45098 
4500 


FÜR ! = S TO R2 STEP 5 4 
HPLOT I X 279 / R2,144 TO I X 279 / R2,149 
NEXT 
HPLOT 2,DF¢@) * .75 + 75 
FOR I = 1 TO F2 
HPLOT TOI x (279 / R2),DF(I) Х „75 + 75 
NEXT 1 

PN = 3 | 
PRINT "DIFFERENCE PLOT " 

PRINT NASE 1) ,NAGS( 25 
PRINT ЕМС1) ,РМС2) 
РЕТМТ 5801),5662) 
END 
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Е 





ӘӘ 
58 
70 


v 

іда 
110 
120 
149 
150 
158 


NEXT 


MAP 0 PROGRAM LISTING 


ШЕ- | : HOME 


РВЇМТ " MAP 8" 
PRINT : PRINT 
PRINT "THIS PROGRAM GENERATES A DATA FILE TO" 
PRINT "BE USED AS INPUT FOR THE TOPOGRAPHIC" 
PRINT "MAPPING SERIES. THE OUTPUT OF THIS" 
PRINT "PROGRAM MAY BE SEEN AS A TARGET" 
PRINT "COVERAGE PLOT. WHEN USED WITH THE" 
PRINT “CONTOUR PLOT, VARIOUS VALUES OF PERCENT" 
PRINT "OF DETECTION MAY BE PLOTTED. 
PRINT PRINT "FOR STANDARD ANALYSES, THE ISOLIN 
ES" 
PRINT "ARE RECOMMENDED TO BE 25, 50 AND 75" 
PRINT "PERCENT." 
INPUT "HIT RETURN TO CONTINUE" ¡As 
HOME 
REM DEFINE CONSTANTS 
ME 1>58б54+С2 = ‚115194 
СЗ - .0й00344:С4 - .819527 
шаг X ATN ( I) :F2 = 50 
0% = CHRS (4) 
MN = 180:MX = ð 
TR = 1389 
REM INPUT VARIABLES 
INPUT "MAME OF PROJECT: ":N$ 
INPUT "SCALE OF PLOT: ";SC 
INPUT "FIGURE OF MERIT :"ıFM 
INPUT "STANDARD DEVIATION OF ESTIMATE: "¡36 
INPUT “HUMBER OF BUOYS IN PATTERN :"sNi 
51 = SC / 40:52 = SC Y 2183 = 2 x SC 
IF S3 < R2 THEN S3 = R2 
REM DIMENSION ARRAYS XXXX 
DIM CEC 46 ,4@) 
DIM PD(S3) ,PQ(S3) 
DIM BxX<N1): DIM BYCN1) 
REM INPUT BUOY POSITIONS 
БШЕ = 1 70 Ni 
БЕТ “INPUT X,Y OF BUOY # "I 
INPUT BXCI) ,BY(D 


I 


(Т 





450 
450 
470 
430 
510 
520 
540 
545 
546 
259 
360 
2/0 
299 


588 


410 
520 
539 


540 


1900 
4799 
5000 
2914 
3012 
5013 
3816 
Su? 
3919 
5020 
5025 
5030 
5449 
5058 
3955 
5060 
5020 
588 
5090 


НОМЕ 

INPUT “PROP LOSS FILENAME {DEL -PLF): ":B$ 
PRINT D$S"ÜPEN ":B$ + "-PLF" 

PRINT D$"READ ":B$ + "-PLF" 

FOR I = @ TO R2 

INPUT PQCI) 

NEXT I 

PRINT D$"CLOSE "ıB$ + "-PLF" 

IF S3 < = 50 THEN 400 

FOR I = R2 + 1 70 2 X SC 
PRTI = 99 

NEXT I 

REM CONVERT PROP LOSS CURVE TO PROBABILITY GF 


DETECTION CURVE 


CEE Р ООЛ ВЕЕ (iy 2X (1 + CI X X t U2 5 
БТ С 532 + Ed X CX 400“ 42 

FOR I = ü TÜ 2 x SC 
ӘБ ЕМ — PULCI) 

ESE ee у ТИЕМ РОК = "EMPL se / 5062: GOTO 
¿58 
ШІ =з 1 = € FN PC == SE Y 28522 

NEXT 1 ' 

с050В 50890 


REM CHANGE FOM OR SIGMA X%xXXx 
INPUT "NEW FOM AND SIGMA Cr/N> 2? 3° 36% 
IF G$ = "N" THEN GOTO 1006 

INPUT "NEW PROJECT NAME: "N$ 

INPUT "MEW FIGURE OF MERIT: “;FM 
INPUT "NEW SIGMA: "356 

GOTO 510 

END 

REM CREATE DATA FILE ххх 

PRINT D$"OPEN";N$ + "-HIRES" 

PRINT D$"WRITE "SNS + "-HIRES" 

PRINT 48 

PRINT 49 

PRINT 488 

PRINT 400 
KT = 21 

СЕЕ ЕМ Ах) = «КТ - 9 % Si 

FOR I - 1 TO 46 


L3 = ЕМ а« 1) 
ӘШЕ = 1 то 90 
і = — FN aco 
KK = 1 


FOR B = 1 ТО М1 
ІСІ > вв: ра — LS = BYCBED 
EIS SUR срез ж LZ + LI. *x Га)» 
ИЕ ИСЕМ ВЕ = CN ~ M:DP = PDM * 15» — PDEM 
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2188 
2118 
5120 
2136 
5140 
5150 
5170 
3250 
ЕИ 
3288 
5268 
5320 
3280 
5464 
9599 


16000 


Ра = 


гр тэ Res OR 
аг PU THEN CECI, J) = PO 
NEXT B 


IEEE CE ID THEN MN = CECI, J? 
IF MX < CELI, J?) THEN MX = CECI, J) 
NEXT 4 
=ÜR J = 1 TO 48 
RINT IMT CIP X CEY(I,J)) Z TR 
NEXT J 
EXT I 
PRINT HN 
PRINT MX 
ШЕЕ CLOSE pit RES" 
RETURN 
PRINT DS"PRHO" 

END 
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MAP I PROGRAM LISTING 


ӘН TEXT НОМЕ 


SES PRINT “ МӨР 1" 


10 
15 
29 
25 
30 
35 
49 
45 


Е 
“2 


55 
59 
28 
25 
т? 
188 
118 
15 


120 


410 


429 
43% 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
НОМЕ 
FEM 


REM 


: PRINT 

"THIS PROGRAM GENERATES A DATA FILE TQ" 
"BE USED AS INPUT FOR THE TOPOGRAPHIC" 
"MAPPING SERIES. THE OUTPUT GF THIS" 
"PROGRAM IS A PROBABILITY OF DETECTION" 
"PLOT. THIS PROGRAM MAY BE USED FOR" 
"LOCALIZATION WITH SINGLE BUDY CONTACT" 
"THE PROGRAM RELIES GN BOTH POSITIVE AND" 
"NEGATIVE INFORMATION. THE FIRST BUDDY" 
"INPUT 15 THE BUOY IN CONTACT." 

: PRINT 

"HIT RETURN TO CONTINUE" ¡As 


DEFINE CONSTANTS 
5854: СЕ (12115174 
.090344:С4 .М1952/ 
45 

4 Х 
CHRS 
180 :MX 
1859 
INPUT VARIABLES 


АТМ 


(4) 
= B 


(1 :R2 = 585 


51 


INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
INPUT 
ТЕ М2 
= 5 
ШЕ S3 
КЕМ 

md Е 
DIMP 


Диме ра 


FOR I 
PRINT 


"NAME OF PROJECT: 
"SCALE OF PLOT: ";SC 
"FIGURE OF MERIT :";FM 
"STANDARD DEVIATION OF ESTIMATE: 
"NUMBER OF BUOYS IN PATTERN :4:М1 
"NUMBER üF BUOYS INB CONTACT : ";N2 
> Nit GOTO 275 

C / 48:82 = SC / 2:5: 
€ R2 THEN $3 = R2 
DIMENSION ARRAYS <**x 

EC 40,48) 

0653) ,РИС53) ,РЕС53) ,РЕС5З) 

DIM BY<N1> 

= 1 TO N2 

"INPUT X,Y OF CONTACT BUOY 


"INS 


Е 


= ž X 


3C 


#";1;: INPUT ^ 


ШЕСІ) (ЕСІ; 


МЕХТ 
IF N2 
REH 

FOR I 
PRINT 


INPUT 
NEXT 


Ni THEN GOTO 454 

INPUT BUOY POSITIONS 

= N2 + 1 TO Mi 

"INPUT X,Y DF BUOY # ";1;" <NO CONTACT) 
BX(1) ,BY{D 

I 
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450 
450 
370 
480 
510 
229 
540 


546 
259 
568 
320 
э?? 


589 


DETECTION CURVE 


GOSUB svu 
GOSUB 5488 
REM CHANGE FÜM CR SIGMA ххх 
INPUT "NEW FOM AND SIGMA 
IF 6% - "N" THEN 0010 {CUM 

INPUT "NENW PROJECT NAME: " NS . 
ПЕРЕН FIGURE OE MERIT: "FM 
INPUT "NEW SIGMA: ";56 
GUTO 514 

END 

REM CREATE DATA FILE X*X* 
EEN EST OPEN SNS + -HLRES"' 
БЕЛИ OE ARRE ORNS + “-HTRES” 
PRINT 4% 

PRINT ай 

PRINT 489 

PRINT 308 
penu 21 

DEF FN ACK) = 
FOUR 12502110 49 
L3 = FN ACI? 
FOR J = 1 TO 46 
L1 = - FN ACJ) 
KK = 1 

FOR B = 1 TO Ni 
L2 = L1 - Bx(B)> :L4 = L3 
IS SOR Ez S Lz - 
ETB < = N2 THEN 


Sim EE 


— ВҮ*В, 
х 14) 
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CY/ND ? 2" 36% 


HOME 

Cee LUSS FILENAME (DEL -PLF): "¡ES 
ЕКШШ EE PEIN PLE” 

ERT UE REGUL SB + "PLP" 

FOR TSU TO R2 

INPUT PHCI? 

NEXT 1 

FERIT ҮГЭЭ SEE ¢ "RLF" 

IF 53 < = 50 THEN ors 

SOI ЕТІ па вт SC 
РӨСІ) = 977 

NEXT I 

REM CONVERT PROPR CüSS CURVE 10 PROBABILITY OF 


+ 
ао? 


ша и = 1 — 20/72 5 £1 + Ll £ x 
шише tes X ХХ 7 22 % 86 x X ^ 45» 

FOR I = ñ TÜ S3 
эж ЕМ = Pec]? 

ШЕ ЕИ = 0 TREN FOCI =: T FEN PSE 7 Sb): 
5548 
БОО — FM PC = "SE 7 50) 

NEXT 1 


га 





5a95 M = INT ¢LN + C5) 

5199 KK = KK Х РЕ(М): GOTO S124 

5119 М = INT (LN + C5) 

5115 КК = KK Х РЕСМ) 

5120 МЕХТ В 

5125 IF MN > KK THEN MN = КК 

5127 IF MX < KK THEN MX = KK 

5130 CECI, J) = KK 

5150 NEXT J 

5170 FOR J = 1 TO 48 

Z25H PRINT INT «ТЕ & CEC1,D> Y TR 

5270 NEXT J 

5290 NEXT I 

5340 PRINT MN 

5370 PRINT MX 

5398 PRINT D$"CLOSE";N$ + "-HIRES" 

5400 RETURN 

3809 Pl = PD(B):P2 = PDC1) 

$010 IF P] = > Р2 ТНЕМ РЕ(Й) + Р.: РЕВ) = 1 - PDC1 
J): GOTO 5039 

$020 PELO) = P2:PFCa) « | - РОСВ) 

5838 FOR I = 1 TO 33 - 1 

5040 Рі = РО<1 ~ D :P2 = POCI):P3 = PDCI + 1) 

$658 IF PI = > P2 AND P1 > P3 THEN РЕСІ) = РІ: GOTO 
5030 

$850 IF P2 = > P1 AND P2 = » P3 THEN PE<I) = P2: GOTO 
5039 

6070 PECI> = PZ 

60308 ІЕ РІ < = P2 AND P1 ¢ = P3 THEN PFCI> = 1 - 
Pi: GOTO $118 

4890 IF P2 < = Рі AND P2 < = P3 THEN PFCID = 1 - 
P2: GOTO 6118 

5100 РЕСІ) = 1 - РЗ 

6110 МЕХТ 

$148 RETURN 

2999 PRINT О®"РЕ#@" 

19994 END 
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285 
286 


POEL 


DEL. 


BLOCK II PROGRAM LINE СНАХСЕ5 


ІМРОТ MN 
INPUT MX 


170,130 


335, 1818 
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CONTOUR II PROGRAM LINE CHANGES 


Балы ПЧРИТ ММ 

256 INPUT Mx 

229  HPLOT 4.9 TO (4.4 X HC - 1»,8 

378  HPLOT TO (4.4 X HC - 12,09 X HR - 5 


ОЕ ISOS 
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DA 


EL | 
: o с И 
эЛОНАН2 8И12 МАЯӘ0ИЯ 11 23901 


L=. E 
= 





— e u 


6,71 s DN RIE DO от 


9 
(I = ан 4 1,116 28 2 
P š 
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